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Abstract. We examine the relevance of several major material-dependent 
parameters to the magnetic softness in iron-base superconductors by first-principles 
electronic structure analysis of their parent compounds. The results are explained 
in the spin-fermion model where localized spins and orbitally degenerate itinerant 
electrons coexist and are coupled by Hund's rule coupling. We found that the difference 
in the strength of the Hund's rule coupling term is the major material-dependent 
microscopic parameter for determining the ground-state spin pattern. The magnetic 
softness in iron-based superconductors is essentially driven by the competition between 
the double-exchange ferromagnetism and the superexchange antiferromagnetism. 
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1. Introduction 



Recently, high transition-temperature (T c ) superconductivity has been observed in a 
number of doped iron-based layer materials P, El |3|, H] near a static antiferromagnetic 
(AF) order [HJ El El [HI E] and with a spin resonance [TQl [HI H2], a pattern exhibited 
previously by the layered cuprate high-T c superconductors. It has been generally 
believed that strong spin fluctuations in two-dimensional (2D) space is at the heart 
of the high-T c mechanism. A proper understanding of the magnetism in parent 
undoped materials thus becomes an essential first step towards understanding the high- 
Tc mechanism. The undoped cuprates are universally described by the 2D Heisenberg 
model. The key character of this model is that a novel resonating- valence-bond spin- 
liquid state competes fiercely with the traditional Neel spin-solid state for being the 
ground state [13], a situation referred to as magnetic softness. Magnetic softness 
becomes even more apparent in the parent compounds of iron-based superconductors 
(FeSCs), since different ground-state AF spin patterns are truly realized — 'collinear' C- 
type in iron pnictides (e.g., LaOFeAs and BaFe2As2 016]) and 'bicollinear' E-type in 
iron chalcogenides (e.g., FeTei-^Se,,; [HE]), as illustrated in Fig. [I] — despite apparent 
similarity in crystal and electronic structures. To elucidate its nature, it is necessary to 
first identify what material-dependent parameter drives the difference in the magnetic 
pattern; this knowledge will put stringent constraints on minimum theoretical modeling 
of FeSCs. 

In literature, the anion height from the iron plane (-2 an io n ) and the ordered magnetic 
moment (m) have been considered as major material-dependent parameters [T6| ITT] , 
since neutron diffraction experiments reported z anion = 1.31, 1.35, and 1.73 A and m = 
0.36, 0.87, and 1.70 hb in LaOFeAs, BaFe2As2, and FeTe, respectively j5j El [TJ IB]. 
First-principles studies qualitatively reproduced these observations [TBI US] and further 
revealed that the C-E magnetic transition can be induced by tuning z an i n in FeTei-^Se^ 
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Figure 1. Schematics of two AF spin orders in the iron plane. Red and blue balls 
stand for spin-up and spin-down Fe atoms, respectively. Dashed lines represent the 
'bicollinear' view. The zigzag horizontal and vertical solid thick lines represent the 
-E-type view [14] following the context of manganites [15] . 
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Figure 2. Comparing the total energies of the C-type and -E-type AF states in 
BaFe2As2 as a function of m obtained by fixed-spin-moment calculations (circles) and 
by varying z an i on (squares). In the latter case, Zanion is marked near the corresponding 
datum. 



[20J. Since varying ,2 an i on also varies m, it is unclear whether it is m or something else 
that determines the C-E transition. Another interesting scenario is that the difference 
in the orbital ordering pattern determines the ground-state magnetic pattern: one of 
the initially degenerate Fe d xz and d yz orbitals gets more populated than the other in 
the C-type pnictides [2D 1221 1231 121] , and Fe d Y z = (d xz - d yz ) / \/2 orbital was stipulated 
to be more populated than d X z = {d xz + d yz )/y/2 in the E-tjpe FeTe [25J. While the 
former was verified by first-principles calculations [22], the latter remains untested. 

The purpose of this paper is to examine the relevance of z anion , m, and orbital 
order to the magnetic softness in the parent compounds of FeSCs by first-principles 
electronic structure analysis. We show that increasing ^ an i on in BaFe2As2 also increases 
m and drives the magnetic transition from C to E, but that varying m and fixing z an ion 
cannot. We further demonstrate that the E type does not exhibit a site orbital order. 
These results are explained by use of a recently proposed spin-fermion model [2] where 
localized spins and orbitally degenerate itinerant electrons coexist and are coupled by 
Hund's rule coupling (2S1 123 I2H1 12H]- We found that the difference in the strength of 
the Hund's rule coupling term is the major material-dependent microscopic parameter 
for determining the ground-state spin pattern in the parent compounds of FeSCs. The 
effort to gain the Hund's rule coupling energy is known to induce the double-exchange 
ferromagnetism [30l l3Tj . Therefore, the magnetic softness in FeSCs is essentially driven 
by the competition between the double-exchange ferromagnetism and the superexchange 
antiferromagnetism in the localized-spins sector. 
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2. First-principles analysis 

We performed first-principles calculations in local spin-density approximation of density 
functional theory with full potential, all-electron, linearized augmented plane wave basis 
implemented in the WIEN2k software package [32] ■ We adopted an eight-Fe-atom unit 
cell for all the calculations and the 10 x 10 x 8 mesh of k points. The energy convergence 
is better than 1 meV per Fe atom. 

In Fig. [2j we compare the total energies of the C-type and E'-type AF states in 
BaFe2As2- First, the effect of -2 anion was investigated. Like in FeTei-^Sea; [20] . increasing 
^anion i n BaFe2As2 also increases m and drives the magnetic transition from C to E 
around -2 anion = 1.6 A. Then, to isolate the effect of m from z an i on , we fixed -2 a nion at the 
experimental position of 1.35 A [2] and varied m using the fixed-spin-moment method. 
We found that varying m alone cannot drive the C-E transition. 

To examine a possible orbital order in the E type, we performed first-principles 
Wannier function analysis [22, 33J on FeTe and got the following Fe density matrix: 
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where the x and y axes point to the nearest-neighbor (NN) Fe atoms. To compare with 
the orbital ordering pattern proposed in Ref. [25], the coordinate system needs to be 
rotated by 45° in the Fe plane. The resulting density matrix is rewritten as 



/ 




3z 2 - r 2 


XY 


YZ 


xz 


X 2 -Y 2 


\ 




3z 2 - r 2 


1.44 


0.00 


0.00 


-0.057 


0.03 






XY 


0.00 


1.31 


-0.057 


0.00 


0.00 






YZ 


0.00 


-0.057 


1.04 


0.00 


0.00 






XZ 


-0.057 


0.00 


0.00 


1.06 


0.071 




\ 


X 2 -Y 2 


0.03 


0.00 


0.00 


0.071 


1.13 


/ 



where the X and Y axes point to the next-nearest-neighbor (NNN) Fe atoms. 
Apparently, orbital polarization within the XZ and YZ orbitals is negligible. Therefore, 
site orbital ordering is not the driving force for the E-type spin order in FeTe. 

Below we explain these results using the spin-fermion model that was shown to be 
capable of providing a unified picture for magnetic correlation and electronic transport 
in the parent compounds of FeSCs [HI [Ml E5]- 

3. The spin-fermion model 

The electrons in the degenerate Fe 3d xz and 3d yz orbitals are treated as itinerant 
electrons, and those in the rest Fe 3c? orbitals as effective localized spins. This leads to 
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Figure 3. Phase diagram of the spin-fermion model for J' = J. Black and red 
arrows correspond to changing z a nion and m in Fig. [2] respectively. The dots marked 
with 1111, 122, and 11 are schematic locations of LaOFeAs, BaFe2As2, and FeTe, 
respectively. 



an effective two-orbital double-exchange model [T4J [26j [27] : 
H = — (ty Cj^^CjY^, + h.c.) 

K ^ - - - 

-yE Qtm^w'Qtm' ' &i + ^2 hiSi ' Sj, (1) 

where C iryfl denotes the annihilation operator of an itinerant electron with spin [i =f or 
I in the 7 = d xz or d yz orbital on site i. €[J 's are the electron hopping parameters, a^i 
is the Pauli matrix and Si is the localized spin whose magnitude is S. K is the effective 
Hund's rule coupling. Jij is the AF superexchange couplings between the localized 
spins; in particular, J and J' are respectively the NN and NNN ones. The filling of the 
itinerant electrons is on average three (one hole) per Fe site, corresponding to the high- 
spin configuration of Fe 3d 6 [22] ■ To the y direction, the d xz -d xz NN hopping integral 
t\\ ~ 0.4 eV and the d yz -d yz NN hopping integral t± ~ 0.13 eV; they are swapped to the 
x direction; by symmetry the NN interorbital hoppings are zero; the NNN intraorbital 
hopping integral t' ~ —0.25 eV for both d xz and d yz orbitals, and the NNN interorbital 
hopping is ±0.07 eV. KS ~ 0.4 — 0.8 eV. JS 2 and J'S 2 are of the same order of 10 
meV. Since the Se anion is located above the center of the Fe plaquette, hybridizations 
via the Fe-As-Fe path give rise to comparable NN and NNN parameters. 

To get a general and simple picture about the magnetic landscaping of the model, 
we compared a variety of static spin orders, such as the ferromagnetic (FM) state and 
the AF states of C-type, E-type, and G-type (i.e., the Neel state where all NN spins are 
antiparallel), with the localized spins treated as Ising spins [HI [36] . Then, Eq. (1) is 
reduced to a system of noninteracting electrons moving in an external potential that is 
—KS/2 and KS/2 at site i when the itinerant electron is spin parallel and antiparallel 
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to Si, respectively. 

The resulting phase diagram of the model is presented in Fig. [3] Changing the spin 
moment S will lead to changing both KS and J'S 2 . This direction of change, as marked 
by the red ray in Fig. |3j will not induce the C-E transition. On the other hand, we 
argued that changing -z an i n corresponds to changing KS, as marked by the black ray 
in Fig. [3] Since the iron atoms communicate with each other through the anions, the 
farther away the anions are, the more isolated the iron atoms are. The isolation of the 
Fe atoms would enhance the local parameter S but suppress the nonlocal parameters 
Jij] therefore, JijS 2 as a whole is much less material dependent than KS. Hence, the 
effective Hund's rule coupling term is decisive in determining the ground-state magnetic 
pattern. (Note that the enhancement of the on-site Coulomb and Hund's rule interaction 
terms by increasing z an i on were also found in the Hubbard model of the 10-fold Fe-3d 
bands for FeSCs [3?].) The results of the spin-fermion model agree well with those from 
first-principles presented in Fig. [2] 

The Hund's rule coupling brings in a blockade mechanism: The electron hopping 
to a neighboring site with opposite spin orientation suffers from the energy barrier 
whose height is of order of KS. Hence, the larger KS is, the more likely the 
neighboring spins are aligned to the same direction. This is the so-called double- 
exchange ferromagnetism [30, 31J, which has been studied intensively and extensively 
in the context of manganites for their colossal magnetoresistance phenomenon. In 
addition, comparable J and J' between the localized spins favor the C type. Hence, 
FeSCs contain intrinsic competition between the double-exchange ferromagnetism and 
the C-type superexchange antiferromagnetism in the localized-spins sector. The E type 
emerges in the intermediate regime as a compromise. 

This model also reproduces the weak orbital polarization within the xz and yz 
(or XZ and YZ) orbitals in the E type because (i) the on-site occupation numbers 
are always the same in both the xz and yz orbitals in the E type from symmetry 
consideration, and (ii) the parameters that can mix these two orbitals to polarize in the 
XZ and YZ representation are negligibly small. 

4. Discussion 

The double-exchange ferromagnetism was employed in the previous proposal that the 
metallic bicollinear antiferromagnetism in FeTe was driven by a site orbital ordering [25J . 
However, the double-exchange effect was treated as being secondary to the Y Z ferro- 
orbital order and the AF spin order along the Y direction; it was used to introduce 
weak FM ordering of those AF Y chains along the X direction. This way actually 
decoupled the whole Fe planar lattice into two interpenetrating sublattices with each 
one exhibiting the C-type AF order on its own. This means that the NN exchange 
coupling strength is weak. On the contrary, in the spin-fermion model Eq. for the 
.E-type FeTe KS ~ 0.8 eV is the leading energy scale; to gain Hund's rule energy is so 
important that the NN Fe-Fe bonding along the zigzag FM chain (see Fig. 1) is strong. 
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This means that if one has to fit the magnetic energy surface of FeTe to the Heisenberg 
model, the resulting NN exchange coupling on average is strong and FM, in agreement 
with neutron scatter data on FeTe [38] . 

The spin-fermion model Eq. ([I]) appears in form similar to that used to describe 
the manganites. These two classes of materials do share some common features, such 
as magnetic softness and large normal-state resistivity. We emphasize two differences 
between FeSCs and manganites. Firstly, the Jahn- Teller distortion energy in FeSCs 
[22J is one order of magnitude smaller than in manganites [33] because the layered 
structure of edge-sharing anion tetrahedra surrounding Fe atoms is much less flexible 
than the network of corner-sharing octahedra surrounding Mn atoms. As a result, unlike 
in the manganites where the orbital degree of freedom is often frozen into an orbital 
order via the cooperative Jahn- Teller effect, there lacks such a locking mechanism in 
FeSCs. A small change in KS could lead to a dramatic change in the orbital order 
status and thus would substantially affect the fit to the Heisenberg model. In fact, 
neutron scattering data on CaFe2As2 [39] and FeTe [38] reported that the C and E 
types were surprisingly well separated in the Heisenberg model parameter space, with 
the leading NN exchange interaction being AF and FM, respectively. This implies that 
while the Heisenberg model is still reasonable for describing magnetic linear response 
near one particular ground state, it cannot capture the essential orbital physics and is 
problematic for describing the general magnetic softness in FeSCs. Secondly, KS ~ 2 
eV in manganites is much larger, driving the system into the FM regime, as shown in 
Fig. 3. With a moderate KS, FeSCs are located near the C-E phase boundary and the 
FM phase becomes an irrelevant high-energy state, providing a necessary environment 
for the formation of singlet superconductivity where the paired electrons have opposite 
spins. 

Regarding superconductivity, this magnetic softness is actually a two-blade sword. 
On the one hand, with its strong ability to undertake electronic reconstruction, the 
system may on the one hand manage to efficiently screen direct electron-electron 
Coulomb interaction and generate appropriate bosonic excitations to glue electron pairs 
in order to form superconductivity. The magnetic softness also implies a large scattering 
rate for electron transport and the observed large resistivity in the normal state. This is a 
favorable feature, according to Homes's law that T c is proportional to both the superfluid 
density and the normal-state resistivity [ID]. On the other hand, a too correlated system 
may find ways to satisfy the competing non-superconducting players first. For example, 
the C X E\_ X region, a mixed C-type and E-type AF, could be formed in accord with 
the change in doping level [TSl E] . Only until none of these players can be adequately 
satisfied could superconductivity show up to relief the high entropy — likely with a low 
superfluid density. While the high-T c mechanism of superconductivity in both FeSCs 
and cuprates remains to be discovered, the present results suggest that FeSCs, though 
closer kin to the manganites than the cuprates in terms of their magnetism, can exhibit 
a quantum phase transition to superconductivity instead of ferromagnetism. 
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5. Summary 

We examine the relevance of the anion height from the iron plan, the ordered 
magnetic moment, and the orbital ordering pattern to the magnetic softness in iron- 
base superconductors by first-principles electronic structure analysis of their parent 
compounds. We conclude that the anion height is an effective tuning parameter and the 
others are not. These results are shown to be compatible with a recently proposed spin- 
fermion model where localized spins and orbitally degenerate itinerant electrons coexist 
and are coupled by Hund's rule coupling. This implies that the difference in the strength 
of the Hund's rule coupling term is the major material-dependent microscopic parameter 
for determining the ground-state spin pattern, and that the magnetic softness in iron- 
based superconductors is essentially driven by the competition between the double- 
exchange ferromagnetism and the superexchange antiferromagnetism. 
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